Bogoliubov-like mode in the Tonks-Girardeau Gas 
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We reformulate ID boson-fermion duality in path-integral terms. The result is a ID counterpart 
of the boson-fermion duality in the 2D Chern-Simons gauge theory. The theory is consistent and 
enables, using standard resummation techniques, to obtain the long-wave-length asymptotics of the 
collective mode in ID boson systems at the T onks-Gira rdeau regime. The collective mode has the 
dispersion of Bogoliubov phonons: u{q) = q\J pU (q)/m, where p is the bosons density and U(q) is 
a Fourier component of the two-body potential. 
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ID boson systems attract much attention 0, 0, 0, 0, 
El 0] in light of recent experiments on cigar-shaped 
atomic traps [l| and gases exposed to linear carbon nan- 
otubes |2j . The history of theoretical studies of ID bosons 
goes back to the celebrated work by Leib and Liniger, 
who found exact integrability of zero-range interacting 
bosons via the Bethe ansatz, for all values of the interac- 
tion strength There are two limiting cases for ID sys- 
tems: the high-density weak-interaction Thomas-Fermi 
(TF) regime Q , where the Bogoliubov energy functional 
[10| and the thermodynamic limit of the Gross-Pitaevskii 
mean-field theory | ll| apply; and the low-density strong 
interaction Tonks-Girardeau (TG) regime of impenetra- 
ble bosons 0G]|. 

Unlike the TF regime, in the TG regime the boson 
wavcfunction is Fermi-like and the fermion-boson duality 
method have been proposed for this regime 0] . The idea 
of this method is to substitute dynamical correlations due 
to strong two-body interactions by statistical correlations 
due to "Fermi" statistics of impenetrable bosons in ID 
0,0- For analytical studies, however, the original Leib- 
Liniger first-quantization approach as well as the boson- 
fermion duality method are complicated since one has 
to solve the Schrodinger equation for the many-particle 
wavefunction, which is a function of a macroscopically 
large number of variables N, the total number of particles 
in the system. Therefore, it is desirable to possess a more 
handy method, e.g., a method based on path-integral 
concepts. 

In this letter we propose such a path-integral approach 
for the exact transformation between ID bosons and ID 
fermions. The method is actually a ID analogue of the 
2D composite particles formalism, which enables map- 
ping fermions on bosons and vise versa in 2D by the cou- 
pling to the Chern-Simons gauge field ^4|. As a result, 
we obtain for the first time the long-wave-length asymp- 
totics of a collective mode in the TG regime. The mode 
turns out to be of a Bogoliubov form. 

We start from a secondary quantized Hamiltonian of a 
homogeneous system of spinless ID bosons, which inter- 



act through a two-body potential U: 
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K + U, 
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K = — I dx^ b {x)(-id x Yi, b {x), 



U 



dxdx'Spb(x)U(x — x)5pb(x) 



where m is the mass; d x — d/dx; Pb(x) — w'tpbix) with 
xpb being the boson operators and 5pb(x) = Pb(x) — pb 
is the density fluctuations where pb is the average boson 
spatial density. 

The many-particle wavefunction in the TF regime ex- 
periences no crucial changes as one particle passes an- 
other, whereas in the TG regime it falls down almost to 
zero as the coordinate of one particle approaches the posi- 
tion of another one. The " fermionized" boson wavefunc- 
tion has zero value if the position of one particle coincides 
with that of another one. Such a reduction of the Hilbert 
space to fermionized wavefunctions is an approximation, 
but is justified far inside the TG regime. 

A fermionized boson wavefunction can be constructed 
from a fermion antisymmetric wavefunction in the fol- 
lowing way 0, El : 



ipb ({xi}) = n si & n fa ~~ d x i}) ■ 



(i) 



In a secondary-quantized language Eq. corresponds 
to the introduction of the new quasi-particle operators 
$t, which are related to by 



4>l{x) = ipf(x) exp ( — in j dx'6(x — x')p(x') 



(2) 



where 9 is the Heaviside unit step function and 

p(x) = ipl(x)ip b (x) = i)\{x)'ij)f{x) 

is the spatial particle density, which has the same form 
in terms of initial bosons and the new quasi-particles. It 
is easy to prove that 

^ f ( Xl )^ f (x 2 ) - e-^ty^Wfixx) 

= \$f( x l)A}(X2)\ = S(xi - X 2 ), 



2 



where A = 0(x% — x 2 ) — 9(x2 — x\) — sign(xi — x 2 ) = 
±1. That is, the new operators satisfy Fermi anti- 
commutation relations so that the quasi-particles are 
fermions. Let us call them composite fermions (CF), af- 
ter their predecessors in 2D. 

To see that the transformation corresponds to 
Eq.JTJ, note that if one starts creating boson wavefunc- 
tion by repeatedly acting on a vacuum state with the 
operators from Eq.(J2|, then the CF operators will pro- 
duce fermionic wavefunction and the exponential phase- 
factors will give the "sign" term in Eq.QJ. 

The kinetic energy and two-body interaction in the CF 
operators' representation take the following forms 



K 



1 

2m 
1 

2~ ; 

with the constraint 



dxipt (x) (— id x + kp + a x ) 2 tpf(x) 



U = — j / dxdx a x (x)U {x — x')a x (x'), 



a x (x) = nSp(x), 



(3) 



where kp = Trp is the Fermi wave- vector. 

In a path-integral representation the constraint J2JI is 
easily incorporated with the aid of a Lagrange multiplier. 



The partition function has the following form so far: 
Z{4>) = J TtyjVtfjDa^S dtc J] 6 ( ^ - p + p) , (4) 



{t,x} 



C 



dx (r f (idt)i>f + p4>) - K(r f ,<Pf) - U(a x ), 



where p = ipjipf an( i the constrained path integration is 
over the statistical field a x and the Grassmann fields ipf 
and which represent the CFs. An external potential 4> 
has also been added to the action. We will use it to probe 
the system, i.e., the density-density correlation function 
is: 



(p(m)p(^,0) = -^)- 1 ^§^ 



(5) 
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The constraint in the partition function (@J can be 
rewritten through the introduction of an auxiliary field 
at as: 

Y[ 8(^- P + p) = j Va t ^ X )e-^ dtd <^'" + p) a \ 

{t,x} 

As a result, the spatial density of the Lagrangian becomes 



(-id x 



2m 



ipf - a t 



(Ox_ 
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2tt 



, f dx'a x (x)U(x — x')a x (x') 



and the integration in Z is assumed now over a time- 
space statistical field a = (at,a x ). 

Before proceeding further, let us outline the connection 
of the proposed boson-fermion transformation to exist- 
ing theories. The transformation @ is reminiscent of 
the inverse boson-fermion transformation in Haldane's 
bosonization approach for ID Fermi liquids |15|. Nev- 
ertheless, Haldane's bosonization is developed for stud- 
ies of low-frequency physics of fermion systems in terms 
of bosons, which, in fact, represent the spatial density 
fluctuations of the fermions. In our case, however, the 
transformation © involves operators of real bosons and 
not density fluctuations. Density fluctuations, instead, 
are represented by the spatial statistical field component 
a x . 

Physically, the proposed model is a ID counterpart of 
Jain's mechanism of attaching flux quanta to 2D parti- 
cles, which leads to the couplin g of the composite objects 
to a Chcrn- Simons gauge field UJ. There are, however, 
several aspects in which the theory proposed differs from 
the Chern-Simons theory, apart from the different dimen- 
sionalities of the systems, (i) The proposed theory is not 



gauge invariant, (ii) Time-reversal and space- reversal 
symmetries are broken separately, though time-space- 
reversibility is present, (m) The coupling to Chern- 
Simons gauge fields in 2D results in an additional mag- 
netic field experienced by the composite objects, whereas 
in our case the coupling to the statistical fields leads to 
a shift of the one-particle kinetic energy dispersion by 
kp in momentum space. In a sense, it looks as though 
the whole system starts moving. This fact is going to 
reveal itself later through a Doppler shift in the response 
function. 

The noninteracting part of the action becomes the sum 
of two Gaussian actions for the CFs and the statistical 
fields governed by the following "bare" CFs' and statis- 
tical fields' propagators respectively (in Fourier space): 



Go\p) = e- 



(p-kp) 2 
2m 

1 

1 v F u(q) 



with u(q) = U(q)/(vpn), where U{q) is the spatial 
Fourier transform of the two-body potential and p = 



3 



a) cr-p 



= 00 + 



b) °^ 




FIG. 1: (a)-(c) Graphical representation of Egs. 16181 . The 
hollow interaction vertices in (c) denote the density compo- 
nents of the two-component density-current interaction ver- 
tex, (d) The excitation spectrum of a ID boson system at 
the TG regime. Besides free particle-hole pair excitations 
there is a collective mode with a Bogoliubov-like dispersion 
uj(q) = qy/pU(q)/m. 



(e,p), q — (to, q). The interaction part of the action is: 



£4 



• + a t )p - a x j, 



where j is the CF current density: 

j = -^-(^/(("^ _ k F)ipf) + ((id x - k F )tp*f)ipf + a x p). 

Due to the coupling to the statistical fields the one- 
particle dispersion is shifted by kp in momentum space. 
One can formally make a substitution p — kf — > p and 
arrive at the ordinary picture of ID fermions at rest. 

As the action is Gaussian in the CF fields, one can 
integrate them out. The integration leads as usual to 
the fermion determinant in the effective statistical fields' 
action. At this point it is important to note that so far 
we made no approximations. Now in order to obtain the 
density-density response of the system it suffices to leave 
in the effective action only the terms quadratic in the 
statistical fields: 



Zctt(<j>) 



Here, the effective action has the following form 

£ cff = a(-q)(D^(q)-fl D (q))a(q) 
-(a + ^(-q)fl p (q)(a + ^(q), 

where H D and H p are the diamagnetic and paramagnetic 
polarization operators respectively. The renormalized 
statistical fields' propagator is given by a Dyson equa- 
tion: 



D(q) = (D-\q)-tl(q)) 



where 



r%) = n P (q) + fl D (q) 



(6) 



(7) 



is the total polarization operator. Using Eq.JSJ one can 
now find the density-density response in the following 
form (see FigQJ: 



(p(-q)p(Q)) = [U P (q)+U P (q)D(q)fl P (q)] t 



(8) 



The advantage of the proposed transformation to CFs 
is that in normal Fermi liquids |l9j we know the exact 
form of the polarization operator in the long-wave-length 
limit |l6( . Due to the so-called loop-cancellation theorem 
(see, e.g., Ref.^3]), in the long- wave-length limit the po- 
larization operator is given by an RPA expression. In 
fact, the loop-cancellation theorem is a consequence of 
gauge invariance, so that one can argue that we can not 
use it in our case. However, far inside the TG regime 
(see below) the gauge invariance is restored so that the 
theorem still applies. In ID one has (q <C k F ) 



il D (q) = 

n p ( Q ) = 





v F 



1 



W 2 - (I — i0+) 2 



— u> 

Dp 



with a/ = uj/(vFq). Substituting these expressions into 
Eq. © one gets the density-density response in the long- 
wave-length limit: 



(p(-q)p(q)) = 



1 



ttof((u/ - I) 2 - («(<?) + 2 - i0+)) ' 



As mentioned, we obtained the Doppler shift since the 
proposed theory lacks time-reversal symmetry. Never- 
theless, as one goes toward the strong TG regime 



u(q) 



00, 



(9) 



the time-reversal symmetry is restored leading to the fol- 
lowing density-density response: 



(p(-q)p{q)) b = 



v F 



(v(q)q - lO+y 



where the momentum dependent " velocity" of the collec- 
tive mode is: 



v(q) = y/pU(q)/m. 



(10) 



The fact that our theory gives reasonable results only 
far inside the TG regime is natural. The fermionization 
of bosonic wavefunctions Q is a good approximation if, 
and only if, the inter-boson repulsion is very strong, i.e., 
if the condition J^J is satisfied. 

The excitation spectrum of ID bosons in TG regime 
(see Fignji) now appears to consist of two kinds of excita- 
tions, just like in Fermi liquids. The first one, originally 
discussed by Lieb 0], is free particle- hole pairs. The 
dispersion of these excitations is given by the density of 
states of the polarization operator n. In the long-wave- 
length limit they are phonons with velocity vp. The sec- 
ond one is the collective mode with velocity JTUJl. This 



4 



collective mode is a surprising result because it has the 
form of Bogoliubov phonons, and in the TG regime the 
Bogoliubov theory is not applicable. 

Another issue is the effective two-body interaction, 
through which the CFs interact. In the limit © the 
gauge invariance is restored and we are left only with 
density-density interactions in the system, because only 
the time-time component of the renormalizcd statistical 
held propagator survives: 

Dtt(q) = nv F u{q)(u;' 2 - \)/{J 2 - u{q)). 

In the static limit (a/ — + 0) the effective interaction 
between CFs is D tt — irvp- This reflects the well- know 
fact that unlike the initial boson system, in the fermion- 
ized case the two-body interaction becomes weak. No- 
tice, however, that in our case with the increase in initial 
inter-boson potential, u(q), the effective two-body poten- 
tial between the CFs does not vanish, as in Ref.|7j, but 
instead reaches a constant value. 

In the limit of high energy transfer, to' 3> u(q), we have 
Dtt(q) = U(q). This is natural, since at high energies the 
screening is relatively weak so that the particles interact 
via the initial potential U(q). 

Finally, in the intermediate region, 1 < u>' < u(q), the 
so-called overscreening effect takes place and the effec- 
tive potential becomes attractive. This fact may have 
interesting consequences. 

In conclusion, we propose a path-integral method for 
statistics transformation in ID, which leads to the cou- 
pling of composite particles to statistical fields, which are 
representatives of density-current interactions. With the 
aid of this transformation we have found a long-wave- 
length limit of a collective mode in TG regime of the 
boson system. The method proposed, however, can be 
used for further study of TG bosons, e.g., the study of 
response of the system near q ~ 2kp , the question of the 
CDW and BCS stabilities, the effect of statistical field 
fluctuations on CF propagator etc. We leave these inves- 
tigations for a later work. 
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